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SUMMARY 

Graphite/epoxy composite thin shell structures were simulated to investigate damage and fracture progression 
due to internal pressure and axial loading. Defective and defect-free structures (thin cylinders) were examined. The 
three laminates examined had fiber orientations of [90/0/±9] s , where 0 is 45°, 60°, and 75°. CODSTRAN, an inte- 
grated computer code that scales up constituent level properties to the structural level and accounts for all possible 
failure modes, was used to simulate composite degradation under loading. Damage initiation, growth, accumulation, 
and propagation to fracture were included in the simulation. Burst pressures for defective and defect-free shells were 
compared to evaluate damage tolerance. 

The results showed that damage initiation began with matrix failure whereas damage and/or fracture progres- 
sion occurred as a result of additional matrix failure and fiber fracture. In both thin cylinder cases examined (defec- 
tive and defect-free), the optimum layup configuration was [90/0/±60] s because it had the best damage tolerance 
with respect to the burst pressure. 


INTRODUCTION 

Graphite/epoxy composite shell structures are used in various structural components such as aircraft fuselages, 
jet engine cowls, pressure vessels, and rocket motor cases. In these applications, it is important to achieve low 
weight, high strength, stiffness, and safety when composite shells are required to withstand significant internal pres- 
sures. This report discusses thin cylindrical shells subjected to internal pressure. Damage initiation, growth, accumu- 
lation, and propagation to fracture were simulated for three cases. 

For a rational design, it is necessary to quantify the defect and damage tolerance of a candidate structure. 
Assessing this tolerance requires the capability to simulate the progressive damage and fracture characteristics of 
composite structures under loading. Critical components of a structure must remain safe but be able to function 
under loading after experiencing some damage. Damage may be caused by an accident, a defect, or an unexpected 
overloading. A structure’s tolerance to damage is quantified by the residual strength (i.e., the additional load- 
carrying ability after damage). 

Another design consideration is the multiplicity of options such as fiber orientation patterns and constituent 
material combinations, choices that result in a great number of parameters that may be varied for an optimal design. 

Design considerations with regard to the progressive fracture of fiber composite structures require a priori 
evaluation of damage initiation and fracture propagation mechanisms under expected loading and service environ- 
ments. Concerns for the safety and survivability of critical components require that the structural fracture resistance 
under loading be quantified. The present approach bypasses traditional fracture mechanics to provide an alternative 
evaluation method, conveying to the design engineer a detailed description of damage initiation, growth, accumula- 
tion, and propagation that would occur in the ultimate fracture of a fiber composite structure. The results show in 
detail the damage progression sequence and structural fracture resistance during different degradation stages. 
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The laminate configuration is a significant design parameter with regard to composite damage tolerance. In 
general, quasi-isotropic laminates have better damage tolerance (ref. 1). However, in many cases when there is no 
damage, a quasi -isotropic laminate may not be the most efficient with regard to structural strength and performance. 
A rational design process requires that the structural damage tolerance for a candidate design be quantified because 
damage initiation and progression characteristics are much more complex for laminated fiber composites than they 
are for homogenous materials. The structural fracture process of a laminated fiber composite depends on many 
parameters: laminate configuration, fiber volume ratio, constituent stiffness, strength, and hygrothermal properties 
and the fabrication process. 

Recent developments in computational simulation have made it possible to evaluate progressive damage and 
fracture in composite structures by assessing the damage initiation and propagation loads. Taken into account 
are the influence of local defects or flaws, through-the-thickness cracks, and the effects of the fabrication process in 
terms of residual stresses. This report demonstrates that computational simulation, with the use of established mate- 
rial modeling and finite element modules, adequately tracks the damage growth and subsequent propagation to frac- 
ture of fiber composite pressurized shells. The present cylindrical shell simulations constitute an extension of the 
work reported in reference 1 in which the defect tolerance characteristics of cross-ply and quasi -isotropic laminates 
were investigated. The current study examines the effects of angle ply orientation on damage tolerance characteris- 
tics and structural fracture modes. The CODSTRAN (Composite Durability STRuctural ANalysis) computer code 
was used to simulate the damage progression in a variety of composite structures: stiffened, adhesively bonded com- 
posite structures (ref. 2), adhesively bonded concentric composite cylinders (ref. 3), damage progression in bolted 
composite structures (ref. 4), progressive damage and fracture of adhesively bonded pipe joints (ref. 5), damage 
tolerance of composite pressurized shells (ref. 6), progressive fracture of fiber composite builtup structures (ref. 7), 
and progressive fracture of composite subjected to simulated Iosipescu shear testing (ref. 8). 

This report describes the application of computational simulation to defect-free and damaged composite cylin- 
drical shell specimens, taking into account damage initiation and/or propagation mechanisms. The purpose was to 
identify the salient material and structural parameters for a design with damage tolerance considerations. 


CODSTRAN METHODOLOGY 

CODSTRAN is an integrated computer code in which three modules are coupled: composite mechanics 
(ICAN), finite element analysis (MHOST), and damage progression modeling. ICAN (Integrated Composite Ana- 
lyzer) is a composite mechanics computer code (ref. 9) that provides the constituent (fiber and matrix) material 
properties from an available data bank and computes the ply properties and the composite properties (effective prop- 
erties) of the laminate in a hygrothermal environment. The code is based on the theory of the micromechanics of 
composites and the classical laminate theory. ICAN can compute ply stresses by using known stress resultants 
(force per laminate thickness, where force can be a bending, a twisting, or a concentrated load). In ICAN, two fail- 
ure criteria were established for the detection of ply failures: (1) the maximum stress criterion in which individual 
ply failure modes are assessed by ICAN using failure criteria associated with the negative and positive limits of the 
six ply stress components in the material directions 1 to 3 (fig. 1): 
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where S^j a represents the ply stress limit (ply strength) in which the ij subscript indicates the stress component and 
the a subscript indicates the sense as tension and/or compression for normal stresses and as ± for shear stress limits 
on the ply; a is the ply stress. The ICAN composite mechanics module computes S^j a stress limits. (2) The modi- 
fied distortion energy (MDE) criterion takes into account combined stresses and is expressed as 
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where a and (3 indicate tensile or compressive stress, S^ lla is the ply longitudinal strength in tension or compres- 
sion, S^ 22 a * s the transverse strength in tension or compression, and 




(l + 4 v n2 - v m) E i22 + v m) E m 

[ E m E m( 2 + v m + v m)( 2 + v m +vm )} /2 


( 8 ) 


The type of failure is assessed by comparing the magnitudes of the squared terms in equation (7). Depending on 
the dominant term in the MDE failure criterion, fiber failure or matrix failure is assigned. If the first squared term in 
equation (7) that corresponds to ply longitudinal tensile or compressive failure is dominant, fiber failure is assigned. 
On the other hand, if one of the other squared terms corresponding to ply transverse tensile or compressive failure or 
to ply shear failure is dominant, matrix failure is assigned. For example, a laminate with configuration [90/0/±75] s 
has ply stresses at the top ply (90°) and is shown in figure 2. In ICAN, the described failure modes of the plies are 
failure due to fiber fracture in tension or in compression; damage due to matrix fracture in tension or in compres- 
sion; and damage due to intralaminar shear fracture. 

MHOST is a finite element computer code (ref. 10) used to solve structural analysis problems. The code can 
perform linear or nonlinear static and dynamic analyses. MHOST has a library containing a variety of elements, 
among which is the four-node shell element used for the present work. (This element is similar to element 75 in the 
MARC structural analysis finite element code.) By supplying the boundary conditions, the type of analysis desired, 
the applied loads, and the laminate properties (using ICAN), MHOST performs the structural analysis. In addition, 
MHOST provides the computed stress resultants to the ICAN code, which then computes the developed ply stresses 
for each ply and checks for ply failure. 

The damage progression module monitors composite degradation for the entire structure and relies on ICAN for 
composite micromechanics, macro mechanics, and laminate analysis. In this module, the overall evaluation of the 
composite structural durability is conducted. 

The CODSTRAN simulation cycle is shown in figure 3. Proceeding clockwise along the left side, one sees the 
constituent material properties (fiber and matrix) provided by ICAN’s data base, the ply properties computed from 
the micromechanics theory, and the laminate properties computed from the laminate theory. These properties in 
conjunction with the finite element mesh, the loads, and the boundary conditions are incorporated in MHOST, which 
performs the structural analysis and provides the computed stress resultants to ICAN (right side of the figure). ICAN 
uses the laminate theory to compute the ply stresses and to check for ply failure. 

The nonlinear structural analysis in the MHOST code is performed in conjunction with an incremental load 
algorithm. The load is increased in small increments (equilibrium positions). In each position, a number of iterations 
(incremental damages) is performed (fig. 4) and the structure is checked for ply failure. If damage is detected, the 
model is automatically updated with a new finite element mesh and new laminate properties; then another finite 
element analysis is performed and the iterations continue until no further damage occurs (equilibrium position). At 
this point, the load is increased and the above procedure repeated until the final failure of the structure. The damage 
progression, fracture, and collapse of the structure are monitored during this procedure. 

The CODSTRAN code is written in FORTRAN 77 computer language for UNIX operating systems at the 
NASA Lewis Research Center. 
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THIN SHELL STRUCTURES 


Damage and progressive fracture were investigated for internally pressurized graphite/epoxy laminated compos- 
ite thin cylindrical shells. The shell structures consisted of eight 0.136-mm (0.00535-in.) plies resulting in a compos- 
ite shell thickness of 1.088 mm (0.0428 in.). Three laminate configurations of [9O/O/±0] S were considered: 0 = 45°, 
60°, and 75°. The 90° plies were in the hoop/circumferential direction of the shells and the 0° plies were in the axial 
direction. Fiber orientations in the ±45°, ±60°, and ±75° (fig. 2) angle plies were given with reference to the axial 
directions. The geometry of the thin cylindrical shells was identical (inside diameter, 305 mm or 12.0 in.; length, 

762 mm or 30 in.). The closed-end cylindrical pressure vessel was simulated by applying a uniformly distributed 
axial tension so that the generalized axial stresses in the shell wall were half those developed in the hoop direction 
(fig. 5). The composite shells were subjected to a monotonically increasing internal hydrostatic pressure until they 
burst. All laminate configurations, defect-free shells and those with through-the-thickness defects, were examined 
(fig. 5). Defects were modeled at a node located at the half-length of the shell. Each finite element model contained 
544 nodes and 512 uniformly sized rectangular elements (fig. 6). Damage initiation, growth, fracture progression, 
and global structural fracture stages were investigated in each case. For both defective and defect-free shells, dam- 
age initiation was by matrix cracking due to transverse tensile stresses in the 0° plies. For the defect-free shells, fiber 
fractures did not occur until just before the structural fracture pressure was reached. The fiber volume ratio was 60 
percent and the void volume ratio was 2 percent of the total volume of the structure. The residual stresses due to the 
processing were taken into account. The cure temperature was 177 °C (350 °F); the pressurization/use temperature 
was 21 °C (70 °F), the moisture content was zero. The fiber and matrix properties were 

obtained from a CODSTRAN data bank of composite constituent material properties. The corresponding graphite 
fibers (AS-4) and epoxy matrix (high-modulus, high-strength (HMHS)) properties are given in tables I and II, 
respectively, and the AS-4/HMHS ply strength is given in table III. Computed results up to global fracture are pre- 
sented for defect-free shells and those with through-the-thickness defects. 


Defect-Free Structures 

Simulations were conducted of defect-free pressurized shells that had three laminate configurations of 
[9O/O/±0] S for which the ply layups were 0 = 45°, 60°, and 75°. The results are summarized as follows: 

The CODSTRAN simulation of [90/0/±45] $ indicated a damage initiation pressure of 0.799 MPa (116 psi). 
Initial damage occurred near the closed ends of the shell and took the form of matrix cracking in selected plies. The 
damage progressed from cracking in the outermost 0° ply to the 45° ply immediately adjacent, to the innermost 
0° ply, and then to the 90° ply on the interior surface. After damage initiation, all plies gradually sustained matrix 
cracking as the pressure was increased to 0.837 MPa (121 psi). After this stage, the pressure was significantly 
increased without activating additional damage modes. The cylindrical shell burst when the 90° hoop plies experi- 
enced fiber fractures at a simulated burst pressure of 4.80 MPa (0.696 ksi). 

The CODSTRAN simulation of the [90/0/±60] s configuration revealed that when the orientations of angle plies 
were changed to ±60°, the damage initiation pressure was increased to 0.965 MPa (140 psi). Similar to the case of 
±45° angle plies, initial damage took the form of matrix cracking in selected plies. At first, the outermost 0° ply 
(ply 2) developed matrix cracking due to ply transverse tensile stresses; the degradation of ply 2 was followed by 
that of ply 3 (the 60° ply immediately adjacent to ply 2). Subsequently, ply 7 (the innermost 0° ply) and ply 8 (the 
90° ply on the interior surface) developed matrix cracks. After damage initiation, all plies gradually sustained matrix 
cracking as the pressure remained virtually constant. After this initial damage growth stage, the pressure increased 
considerably without activating additional significant damage modes. The cylindrical shell burst when ply 1 (the 
outermost 90° hoop ply) experienced fiber fractures that resulted in a rapid propagation of damage to structural frac- 
ture. The simulated burst pressure was 9.42 MPa (1.366 ksi). 

The computational simulation of the [90/0/±75] s configuration indicated a damage initiation pressure of 
1.062 MPa (154 psi), which was the highest pressure of the three defect-free shells investigated. Damage initiation 
was by matrix cracking due to ply transverse tensile stresses in ply 2 (the outermost 0° ply). Immediately after 
matrix cracking occurred in ply 2, plies 3,7, and 8 experienced degradation. Failure modes were in-plane shear in 
ply 3, transverse tensile in ply 7, and transverse tensile and in-plane shear in ply 8. Damage growth in the other plies 
occurred as the pressure was kept virtually constant. After the initial matrix degradation stage, the pressure was 
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increased without activating additional damage modes. The cylindrical shell burst when ply 2 (the outermost 0° axial 
ply) experienced fiber fractures that caused rapid damage propagation to the other plies and then structural fracture. 
The simulated burst pressure was 7.83 MPa (1.136 ksi). It is noteworthy that structural fracture was initiated by the 
tensile failure of axial fibers in this case whereas in the previous two defect-free shells, structural fracture com- 
menced by the fracturing of the hoop plies. 

In general, overall structural damage may include individual ply damage and through-the-thickness fracture of 
the composite laminate. CODSTRAN is able to simulate varied and complex composite damage mechanisms by an 
evaluation of the individual ply failure modes and the associated degradation of laminate properties. The type of 
damage growth and the sequence of damage progression depend on the composite structure, loading, material prop- 
erties, and hygrothermal conditions. A scalar damage variable, derived from the total volume of the composite mate- 
rial affected by the various damage mechanisms, is also evaluated as an indicator of the level of overall damage 
induced by loading. This scalar damage variable is useful for assessing the overall degradation of a given structure 
under a prescribed loading condition. The rate of increase in the overall damage during composite degradation may 
be used as a measure of the structural propensity for fracture. Figure 7 presents a comparison of the overall damage 
progression histories for the three defect-free cylindrical shells. Damage initiation pressures and progression patterns 
are similar. However, burst pressures are strongly influenced by the laminate configuration. The [90/0/±60] s defect- 
free laminate had the highest burst pressure. On the other hand, the [90/0/±75] s defect-free laminate had the highest 
damage initiation pressure. The [90/0/±45] s defect-free laminate had the lowest burst pressure as well as the lowest 
damage initiation pressure. 


Defective Structures 

Three additional simulations were conducted on shells that were modeled to have a through-the-thickness defect 
in the form of an existing 95-mm- (3.75-in.-) long thin axial slit superimposed on a 60-mm- (2.36-in.-) long circum- 
ferential slit at the half-length of the shell (fig. 6). Including a circumferential as well as an axial slit allowed self- 
similar crack propagation in the axial and circumferential directions. A summary of the computational simulation 
results for the through-the-thickness defective thin shells follows: 

For the [90/0/±45] s quasi -isotropic shell, damage initiation began by matrix cracking adjacent to the defect as 
soon as an initial pressure of 0.172 MPa (25 psi) was applied. However, damage growth remained localized adjacent 
to the defect until a 0.745-MPa (108-psi) pressure was reached. After damage initiation, all plies gradually sustained 
matrix cracking as the pressure was increased to 0.837 MPa (121 psi). The damage accumulation stage (by matrix 
cracking) was similar to that of the defect-free shell. However, after the matrix cracking stage was completed, the 
damage propagation to ultimate fracture was much more rapid than that of the defect-free shell. Damage propagation 
by fiber fractures was concentrated near the defect. Through-the-thickness structural fracture progression occurred at 
1.58 MPa (0.230 ksi), causing the cylindrical shell to burst. For the [90/0/±45] s laminate, the axial crack orientation 
was the more critical. Figure 8 shows the finite element model as the simulated defective shell burst. 

Similar to the ±45 quasi-isotropic defective shell, damage initiation in the [90/0/±60] s shell began by matrix 
cracking adjacent to the defect as soon as the initial pressure of 0.172 MPa (25 psi) was applied. Initially, matrix 
failures at the circumferential slit tips commenced by transverse tensile failures in ply 8, and matrix failures in the 
axial slit tips commenced by transverse tensile failures in ply 7. As in the 45° angle plies, damage growth remained 
localized adjacent to the defect until a 0.91 1-MPa (132-psi) pressure was reached. Substantial matrix degradations 
away from the defect also occurred as the pressure was gradually increased from 0.91 1 to 1.020 MPa (132 to 
148 psi). Fracture progression began at 1.765 MPa (256 psi) by the extension of the axial slit. As the pressure in- 
creased, the circumferential slit also became extended. Structural fracture occurred at 1.90 MPa (0.275 ksi), causing 
the cylindrical shell to burst. 

Damage initiation in the [90/0/±75] s shell began by matrix cracking adjacent to the defect as soon as the initial 
pressure of 0.172 MPa (25 psi) was applied. Matrix failures at the circumferential slit tips were initiated by trans- 
verse tensile failures in ply 8. As in the 45° angle plies, damage growth remained localized adjacent to the defect. 
The ultimate structural fracture stage began as the pressure was increased to 0.621 MPa (90 psi). Fracture progres- 
sion was mainly by extension of the circumferential slit. Structural fracture occurred at 0.72 MPa (0.104 ksi), caus- 
ing the cylindrical shell to burst by a circumferential fracture. 
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Figure 9 compares the overall damage progression histories for the three defective cylindrical shells. Damage 
initiation pressures and progression patterns are similar. However, burst pressures and structural fracture modes are 
strongly influenced by the laminate configuration. The [90/0/±60] s laminate has the highest burst pressure, as is the 
case for the defect-free laminates. On the other hand, the (90/0/±75] s defective laminate has the lowest damage ini- 
tiation pressure, unlike the defect-free shells. The [90/0/±45] s defective laminate has intermediate burst and damage 
progression pressures. An examination of structural fracture modes indicates that for the [90/0/±451 s laminate, burst- 
ing occurs at the commencement of fracture propagation along the axial direction. For the [90/0/±60] s defective 
laminate, fracture also commences along the axial direction; however, a circumferential extension of the defect 
occurs at the same pressure level. On the other hand, for the [90/0/±75] s defective shell, fracture initiation and 
progression occur mainly in the circumferential direction. The circumferential fracture mode of the [90/0/±75] s 
laminate is considerably less damage tolerant. Table IV summarizes the burst pressures for the defect-free and 
defective shells. If damage tolerance is defined as the ratio of the defective to the defect-free burst pressures for each 
laminate, the [90/0/±45] s laminate would have a damage tolerance ratio of 33 percent; the [90/0/±60] s laminate, a 
ratio of 20 percent; and the [90/0/±75} s laminate, a ratio of 9.2 percent for the assumed initial damage. However, if 
all burst pressures had been normalized with respect to the burst pressure of the defect-free [90/0/±60] s laminate, the 
damage tolerance ratio of the [90/0/±45] s laminate would drop to 17 percent and that of the [90/0/±75] s laminate 
would be reduced to 7.6 percent. These results indicate that the [90/0/±60] s laminate has the best damage tolerance 
features in this application because of its well-balanced fiber reinforcement pattern. 

Another perspective on damage progression characteristics may be obtained by quantifying the structural resis- 
tance to damage propagation. The global damage energy release rate (DERR) is defined as the rate of work done by 
external forces during structural degradation with respect to the damage produced. The DERR can be used to evalu- 
ate structural resistance to damage propagation at different stages of loading. Figure 10 shows the DERR as a func- 
tion of the applied pressure for the defective [90/0/±60] s specimen. The DERR for damage initiation is relatively 
small, indicating a low resistance to damage initiation by matrix cracking. However, after the damage initiation 
stage, the DERR steadily increases, which indicates a greater structural resistance to damage propagation during a 
damage stabilization stage prior to global fracture. The maximum DERR level corresponds to the completion of the 
matrix cracking stage prior to the occurrence of ply fiber fractures. 

Figure 1 1 shows the DERR levels for the defective [90/0/±75] s shell, for which the very low levels during the 
initial stages of damage progression indicate the lack of structural resistance to damage propagation and predict the 
sensitivity of damage progression to variations in the constituent material properties. 


CONCLUSIONS 

The computer code CODSTRAN (Composite Durability STRuctural ANalysis) was used to evaluate the dam- 
age growth and propagation to fracture of cylindrical composite shells. The conclusions follow: 

1 . The CODSTRAN simulation demonstrates the progression of damage growth and subsequent propagation to 
fracture for defect-free and through-the-thickness defective composite shells subjected to internal pressure. 

2. Through-the-thickness defects have a significant effect on the structural fracture pressure for internally pres- 
surized thin shells. However, the damage initiation stage begins at low load levels for all simulated laminate con- 
figurations. 

3. The damage stability after a major accident can be evaluated to serve as an index of the structural damage 
tolerance. 

4. Computational simulation, with the use of established composite mechanics and finite element modules, can 
be utilized to predict the influence of existing defects, loading, constituent properties, and residual stresses on the 
safety and durability of composite structures. 

5. The demonstrated procedure is flexible and applicable to all types of constituent materials, structural geom- 
etry, and loading. Hybrid composites and homogeneous materials as well as binary composites can be simulated. 

6. The reduction in the ultimate pressure and the variation of the structural fracture mode are predictable with 
reference to the location and the size of the defect or the damage. 

7. The procedure used herein demonstrated that fracture toughness parameters such as the structural fracture 
pressure and damage progression characteristics are identifiable for fiber composite structures with defects. 
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TABLE I. — AS -4 GRAPHITE FIBER PROPERTIES 


Number per end 

Diameter, mm (in.) 

Density, kg/m 3 (lb/in. 3 ) 

Normal modulus, GPa (psi) 

Longitudinal 

Transverse 

Poisson's ratio 

V12 

V 23 

Shear modulus, GPa (psi) 



Thermal expansion coefficient, /°C (/°F) 

Longitudinal 

Transverse .. 

Heat conductivity, J-m/hr/m 2 /°C (BTU-in./hr/in. 2 /°F) 

Longitudinal 

Transverse 

Heat capacity, J/Kg/°C (BTU/lb/°F) 

Strength, MPa (ksi) 

Tensile 

Compressive 


10 000 

0.00762 (0.300x10 3 ) 
4.04 xlO -7 (0.063) 

227 (32.90X10 6 ) 

13.7 (1.99x10*) 

0.20 

0.25 

13.8 (2.00x10^) 

6.90 (l.OOx 10 6 ) 

1.0 xltr 6 (-0.55xl0 6 ) 
1.0x10* (-0.56x1 O' 6 ) 

43.4 (580) 

4.34 (58) 

712 (0.17) 

3723 (540) 

3351 (486) 


TABLE II. — HMHS EPOXY MATRIX PROPERTIES 


Density, kg/m 3 (lb/in. 3 ) 3.40x1 0" 7 (0.0457) 

Normal modulus, GPa (ksi) 4.27 (629) 

Poisson's ratio, v 0.34 

Coefficient of thermal expansion, /°C (/°F) 0.72 (0.4 xlO"*) 

Heat conductivity, J-m/hr/m 2 /°C (Btu-in./hr/in. 2 / 0 F) 930 (1.25) 

Heat capacity, J/kg/°C (Btu/lb/°F) 738 (0.25) 

Strength, MPa (ksi) 

Tensile 84.8 (12.3) 

Compressive 423 (61.3) 

Shear 148 (21.4) 

Allowable strain 

Tensile 0.02 

Compressive 0.05 

Shear strain 0.04 

Torsional... 0.04 

Void conductivity, Lm/hr/m 2 /°C (Btu-in./hr/in. 2 /°F) 16.8 (0.225) 

Glass transition temperature, °C (°F) 216(420) 


TABLE III.— PLY STRENGTH OF AS^t/HMHS 
[Fiber direction, parallel to ply 1 material axis; 
tension, T; compression, C; material axes of 
ply, 1 to 3.] 


Ply stress component 

Strength 


MPa 

ksi 

S*nx 

1930.30 

280 

o 

1475.85 

210 

S{22T 

91.38 

13 

S<22C 

228.27 

33 

s n2 

65.57 

9.5 

S<23 

59.98 

8.7 


TABLE IV — BURST PRESSURES 
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Laminate configuration [90/0/±75] s 







Displacement 


Figure 4. — CODSTRAN load incrementation. 
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Boundary conditions 


AT X = 0, U x = U y = U z = 0 
AT X = L, U y = U z = 0 


274 



272 

Geometry of crack 


Figure 5. — Applied loads and boundary conditions of thin shell structures. 
Cylinder length, L, 762 mm (30 in.); inner radius, r, 152.4 mm (6 in.); lam- 
inate thickness, t, 1 .088 mm (0.0428 in.); graphite/epoxy fiber volume 
fraction, 60 percent. 
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Damage 



Figure 6. — Finite element mesh of graphite/epoxy composite thin shell. 



Normalized pressure (max. pressure, 9.42 MPa (1.37 ksi)) 

Figure 7. — Damage progression in defect-free graphite/epoxy thin cylindrical shells [9O/O/±0], 


2 





Figure 8. — Burst of defective shell [90/0/±45] s . 



Figure 9. — Damage progression in defective graphite/epoxy thin cylindrical shells [9O/O/±0] S . 
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Normalized internal pre< 


Figure 1 0. — Damage energy release rate (DEI 
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